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We derive Casimir-Polder and van der Waals potentials of one or two atoms with diamagnetic 
properties in an arbitrary environment of magnetoelectric bodies. The calculations are based on 
macroscopic quantum electrodynamics and leading-order perturbation theory. For the examples of 
an atom and a perfect mirror and two atoms in free space we show that diamagnetic dispersion 
potentials have the same sign as their electric counterparts, but can exhibit quite different distance 
dependences. 



PACS numbers: 34.35.+a, 34.20.-b 12.20.-m, 42.50.Nn 



I. INTRODUCTION 

The issue of electrodynamics in magnetodielectrics is 
a long-standing one These range from the form and 
propagation of Maxwell fields in such media, to their ab- 
sorption and emission characteristics, as well as to the 
numerous other optical properties that they exhibit. One 
curious aspect, recently confirmed by fabrication of peri- 
odic arrays of thin metallic wires @t3; is the existence 
of left-handed materials, so-called because the electric, 
magnetic and direction of propagation vectors form a 
left-handed triad, resulting in a metamaterial with neg- 
ative refractive index. Another interesting aspect is the 
nature of inter-particle interactions in media with elec- 
trical permittivity, magnetic permeability, or both. This 
has its origins in the work of Casimir [5|, subsequently 
spawning a wide variety of related phenomena catego- 
rized under Casimir-like effects [u-9], and including the 
Casimir-Polder (CP) dispersion potential [Icj . 

The understanding and computation of van der Waals 
forces (vdW) for chemical, physical and biological sys- 
tems is of fundamental importance and widespread in- 
terest as they manifest in vacuum, gaseous or condensed 
phases. These vary from semi-classical treatments yield- 
ing the London dispersion formula, to quantum field the- 
oretic approaches where radiative effects are properly ac- 
counted for and result in Casimir energy shifts. Perhaps 
the best-known example of a fundamental theory that 
automatically includes retardation is quantum electrody- 
namics (QED). Two often- used versions include macro- 
scopic [lll - ITij and microscopic [H, EH QED, their pre- 
fixes aptly describing their general range of applicability. 



In the former variant, the body- assisted electromagnetic 
field is evaluated by quantizing the radiation field and 
the dispersive and absorptive medium. This enables the 
body-induced atomic energy shift to be computed, which 
is interpreted as the potential of the force. From this, 
pair and many-body interaction energies can also be cal- 
culated. This has been successfully carried out for both 
electrically polarizable [ItJ and paramagnetically suscep- 
tible systems |l8[ . Effects of an intermediate medium can 
also be accounted for [TH, H(| • 

Microscopic Coulomb gauge QED constructed via 
quantization of the free electromagnetic field has also 
been used successfully to evaluate atom-field and atom- 
atom interactions in a vacuum. In the latter situation 
this has included re-calculation of the CP potential, and 
its extension to higher multipole moment contributions, 
such as magnetic di pole, e lectric quadrupole and dia- 
magnetic couplings |2l| - |27t . Unexpected results ensue 
such as the discriminatory nature of the interaction be- 
tween two chiral (optically active) molecules; the repul- 
sive form of the ground state dispersion energy shift 
between an electrically polarizable atom and a param- 
agnetically susceptible one; and that in the near-zone, 
the electric-diamagnetic contribution to this potential 
is larger than the corresponding limit arising from the 
electric-paramagnetic term (23j. 

In this paper a systematic study is performed on the 
CP potentials of a diamagnetic atom and its vdW po- 
tentials with another, and with an atom that is either 
electrically or paramagnetically polarizable in a medium 
comprised of magnetodielectric bodies using macroscopic 
QED theory. Interest in contributions arising from dia- 
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magnetic coupling has been due to their importance 
when computing hi ghly accurate potentials for alkali 
metal atom dimers |28|], where it has been found that 
the electric-diamagnetic and paramagnetic-diamagnetic 
terms can be larger than the electric-paramagnetic and 
paramagnetic-paramagnetic energy shifts. The article 
is organized as follows. Section [TT] briefly describes the 
quantized body-assisted Maxwell field operators in the 
medium, their expression in terms of the Green's- function 
solution of the Helmholtz equation, and the writing of 
the total system Hamiltonian. Electric and total mag- 
netic contributions to the single atom CP shift, and for 
an atom placed in front of a perfectly reflecting mirror 
are obtained in Sec. IIIII Presented in Sec. |IV] are dia- 
grammatic perturbation theory results for the dispersion 
pair potential when at least one of the atoms is diamag- 
netic. Corresponding free-field interaction energies are 
also given, and these are briefly compared with results 
obtained via microscopic QED. In all of the cases ex- 
amined, the asymptotic limits of the energy shift in the 
near- and far-zones corresponding to non-retarded and 
retarded regimes are found, and compared with previ- 
ously obtained limiting forms for electric, paramagnetic, 
and electric-paramagnetic interactions. A short sum- 
mary is given in Sec. |V] 



II. QUANTISATION SCHEME 

Consider one or two atoms (or molecules) with internal 
Hamiltonians (£ = A, B) 



(1) 



Within the multipolar coupling scheme, the interaction 
of each atom with the body-assisted electromagnetic field 
is given by [29| 



(6) 



Here the three terms are the electric, paramagnetic and 
diamagnetic interactions, 



H, 



H?j? = -rh ( B(r^ 



Tjd 

AfF — 



\ - c 



« F ^ 8m a 



[r a xB(r 6 )Y 



(7) 
(8) 

(9) 



with £i£ and rh^ being the respective atomic electric and 
magnetic dipole operators and q a , m a and f a denoting 
the charges, masses and positions relative to the center of 
mass of the particles contained in the atoms. Note that 
electric quadrupole and even octupole contributions can 
easily be included in the formalism; they have recently 
been shown to affect the CP potential of Rydberg atoms 
close to surfaces (30j . 

Introducing the atomic diamagnetisability operator as 



E 



4m Q 



(rll-r a r a ), 



(10) 



and using the identity [a x b} 2 = b ■ (a 2 \ — aa) ■ b, the 
diamagnetic interaction Hamiltonian may be cast in the 
form 



(11) 



The ground-state diamagnetisability of an atom is given 
by the expectation value 



which are placed at positions within an arbitrary ar- 
rangement of linearly, locally and isotropically respond- 
ing magnetoelectric bodies, described by a Kramers- 
Kronig consistent permittivity e{r,uj) and permeability 
/i(r, w). Upon introducing bosonic variables f\(r, and 
f\{r,uj), which are creation and annihilation operators 
for the elementary electric (A = e) and magnetic (A = m) 
excitations of the system of bodies and electromagnetic 
field and obey the bosonic commutation relations 



[f x (r,u)Jl(r',u / )}=6(r-r')6 X y6(u J -u / ) (2) 



and 



[A(r, w ),/ V (r', w ')] = [fl(r,u)Jl(r',cj')] = 0, (3) 
the body-field Hamiltonian takes the form [TJj 



Hp 



E 



d d r 



durkjfl(r,oj)-fx(r,u). (4) 



The ground-state of Hp can obviously be defined by 

fx(r,u)\{0}) = VA,r,«. (5) 



(0£|f£l-f a F o |0*) 



V ~ ' 17 ' 'P a )\ = Pp, 



(12) 



where the second line holds for isotropic atoms. 

The total Hamiltonian of the atom(s) interacting with 
the electromagnetic field in the presence of the bodies 
takes the form 



H 



E ^ 

t=A,B 



JL +H F + H iF . (13) 

t,=A,B 



The electric and magnetic field operators can be ex- 
panded in terms of the bosonic operators ft(r, oS) and 
f\(r,w) as 



E{r) = Y I d '' r ' 

A— e,m 

+ H.c 



B(r) = E / dV 
+ H.c. 



du)G x (r,r',u])-fx(r',u>) 



(14) 

VxG A (r,r»./ A (r', W ) 
(15) 
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The expansion coefficients G\ are related to the classical 
Green tensor G by 



G e (r,r',oj) = Ime(r', oj) G(r,r',oj), (16) 



ire 



~ i i \ .oo I h lmu(r', oj) r „, , ,,t 

G m r, r', oj = i - J | [) \' V'x G(r , r, w) ; 

c y 7re |/4 r 

(17) 

and the Green tensor is the unique solution to the 
Helmholtz equation 



„ 1 _ oj , . 

V x — V x - — e(r, oj) 

fi(r,uj) cr 



G(r,r',uj) = S(r — r') 
(18) 



together with the boundary condition 
G(r,r',w) — > for \r — r'\ 



(19) 



Just like the permittivity and permeability, the Green 
tensor is an analytic function in the upper half of the 
complex frequency plane and fulfils the Schwarz reflection 
principle 

G(r,r',-w*) = G*(r,r',oj). (20) 
In addition, it obeys the Onsager-Lorentz reciprocity 

G(r',r,w) = G T (r,r',w) (21) 
and the useful integral relation 



]T yd 3 S G A (r, S , W )-Gf(r', S , W ) 



X—e.m 



^ W 2 ImG(r,r» (22) 

7T 



holds. 



III. CASIMIR-POLDER POTENTIAL OF A 
SINGLE ATOM 

In this section, we calculate the CP force on a single 
atom in the presence of magnetoelectric bodies, discard- 
ing the label A wherever possible. According to Casimir 
and Polder, the force can be derived from the associated 
CP potential U{r A ) 



-VU(r A ), 



(23) 



which in turn can be identified as the position-dependent 
part of the energy shift AE due to the atom-field cou- 
pling 



A. Perturbation theory 

We assume the atom-field system to be prepared in the 
uncoupled ground state |0 J 4)|{0}) and calculate the en- 
ergy shift due to the atom-field coupling within leading- 
order perturbation theory. Let us first study a purely 
diamagnetic atom, whose CP potential is due to the first- 
order energy shift associated with the diamagnetic part 
of the atom-field interaction (fTTj) . 

AE = ({0}\(0 A \ - \ B{r A )-p d -B(r A )\b A )m). (25) 

Normally, CP energy shifts only arise in second order 
perturbation theory. However, due to the diamagnetic 
interaction Hamiltonian being quadratic in the mag- 
netic induction field, already the first perturbation or- 
der contributes. It can be easily evaluated by using the 
magnetic- field expression (|15p . the bosonic commutation 
relations ([2]) and ((SJ, and the integral relation (|22|) . re- 
sulting in 

AE= f^o f dwtr [^. ImGmm(r ^ rA;W )] (26) 
Z7T Jo 



with 



G mm (r, r', uj) = Vx G(r, r\ oj) X (27) 



For an atom in a free-space region, the CP potential 
can be extracted from this energy shift by separating the 
Green tensor into its bulk and scattering parts, 

G(r,r',oj) = G (0) (r,r',cj) + G (1) (r,r',cj), (28) 

and discarding the constant energy shift associated with 
the translationally invariant bulk part by making the re- 
placement G i-> G^. The result can be simplified by 
converting the integral over real frequencies to another 
along the imaginary axis in the complex frequency plane. 
To this end, we first write 

/•OO POO 

/ dwImG(r,r',w)=Im / dojG(r , r' , oj) . (29) 
./o Jo 

The integral on the right hand side can be replaced by 
an integral along the positive imaginary axis (oj > i£) 
plus a vanishing integral along infinite quarter-circles via 
Cauchy's theorem. Using the fact that G(i£) is real- 
valued for a real £, as can be inferred from Schwarz re- 
flection principle (f20)). one finds 



Ud (r A ) Jj£ /°detr[/3 d .G« n (r A ,r A ,iC)] 

M /ll^trG^A.r^). (30) 
27r Jo 



U(r A ) = AE(r A ). 



(24) 



The second line in Eq. (|30| holds for isotropic atoms. 
For an atom with an additional nontrivial electric and 
paramagnetic response, the respective electric and para- 
magnetic interactions in the coupling Hamiltonian © 
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also need to be taken into account. As shown previously, 
they give rise to second-order energy shifts such that the 
electric and paramagnetic CP potentials are given by [Hi ] 



pi r°° 

U e (r A ) = — detr[a(iO-G«(r A ,r A ,iO] 

^7T£0 JO 

= / d£a(i£)trG«(r A ,r- A ,i£), (31) 

27T£ Jo 

U P (r A ) = ^ /°dOr[/3 p (iO-G^(^,r A ,iO] 
Z7r Jo 

= ^° r d ^f(iOtrG« m (r A ,r A ,iO (32) 
2tt Jo 



with 



G ee (r,r,w) = — 2-G(r,r,w) 



(33) 



and 



a(uj) = lim — - — t- 



2 ^ ^A^A VA 



o ft {u! a ) 2 — lj 2 —iue 
2 ^ 



= lim — — r 

£ ^o 3ft ^ uj\ 



(3 p (uj) = lim-V 

£->0 ft ^ 



(o;^) 2 — w 2 — iuje 
u> A m A k m A ° 



art —ui^ — xuje 



I = a(u)l, (34) 



lim — - — r 

e^o 3ft ^ uj\ 



u)" A \m 



Ofc 1 2 



3ft ^ (w A ) 2 — uj 2 — iuie 



1=I3 p (uj)1 (35) 



= {E\-E%)/K, n°2 



Ok _ 



•\P-A\k), rri 



Ok _ 



\m A \ 



denoting the polarisability and paramagnetisability of 
the atom, respectively. Introducing the total magneti- 
sability 



= [^H+/3 d ]l=/3(c)l, 
the magnetic part of the CP potential reads 
U m (r A ) = U p (r A ) + U d (r A ) 

ft/LtQ 



(36) 



2tt 
ft/^o 



d£ tr[p(iO-G$ m (r A ,r A ,iO] 

O 

d^(iOtr[G^(r A) r^,iC)] (37) 



2vr jo 

and the total CP potential is given by 

U(r A ) = U e (r A ) + U m (r A ) 



(38) 



We have thus generalised previous results for the CP 
potential of an atom with electric and paramagnetic 
properties to one that also exhibits nontrivial diamag- 
netic properties. It is found that despite the different 
interaction terms and perturbative orders (first order in- 
stead of second order), the extension to a diamagnetic 
atom can be obtained formally by including the dia- 
magnetic contribution in the magnetisability, f3 p (u) t— > 



/3(w) = /3 p (ui) + (3 d . In particular, the local-field cor- 
rected potentials for atoms embedded in a medium as 
derived in Ref. [l8[ remain valid with this replacement. 

By introducing a e (cj) = a(ui), a m (u) — f3(ui)/c 2 , the 
electric and magnetic parts of the CP potential can be 
given in the compact notation 



U x {r A ) 



2vre 



2tt£ 



/ d£te[a A (iO.G&(r A) r >1 ,iO] 
Jo 

poo 

/ d£a A (i£)tr G<&(r A ,r A ,%) (39) 
Jo 



(A = e, m). 

There are two important differences between the dia- 
magnetic and the paramagnetic magnetisabilities which 
will have an impact on the associated potentials. Firstly, 
the diamagnetisability has an opposite sign with respect 
to the paramagnetisability, which is a consequence of the 
Lenz rule. Secondly, in contrast to the paramagnetisabil- 
ity which obeys the usual Kramers-Kronig relations, the 
diamagnetisability is independent of frequency. 



B. Application: Atom in front of a perfectly 
reflecting mirror 

Let us consider an isotropic atom at distance z A from 
a perfectly reflecting planar mirror. The magnetoelec- 
tric properties of the mirror are characterised by e = oo 
(/i = oo) for a perfectly conducting (infinitely permeable) 
plate. The Green tensor reads [31j 

G«(r,r',i£) 

= ± ^j ^-^-^H^-^) (40) 



(q± = q± ibe z , q±e z , q = \q\, b = y/q 2 + £ 2 /c 2 ) with 
the upper (lower) sign corresponding to a perfectly con- 
ducting (infinitely permeable) plate and the polarisation 
vectors e s and e p being defined by (e q = q/q) 



(-ige z =p be g)- 



(41) 



Evaluating the double curl of the Green tensor, we find 
that G^r.r'.if) is equal to-£ 2 c - 2 G (1) (r,r',i£). Sub- 
stituting this into Eq. (|30|) and carrying out the q- 
integral, one finds 



U d (z A ) = 



l6ir 2 z A 



d£e 



-2za£/c 



(42) 



After performing the ^-integral, we find an attractive (re- 
pulsive) CP potential 



U d (z A ) 



± 



3ft^ c/3 d 
32ttz^ 



3ft/^oc i 
32^iT ^ 6r 

71 a£A 



(43) 
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of a diamagnetic atom in front of a perfectly conducting 
(permeable) plate. It is given by a single l/z\ power law. 

In Table U we compare this result with the known find- 
ings for electric and paramagnetic atoms. We recall that 
electric and paramagnetic atoms interact with conduct- 
ing and permeable plates according to an 'equals-attract, 
opposites-repel' rule: An electric plate attracts electric 
atoms while repelling (para)magnetic atoms, with corre- 
sponding results for a (para)magnetic plate. In contrast 
to this, diamagnetic potentials carry a sign that is op- 
posite to their paramagnetic counterparts. This is due 
to the Lenz rule as encoded in the minus sign in the 
diamagnetic magnetisability (fi"2j) . The diamagnetic CP 
potential thus has the same sign as the corresponding 
electric potential. 

Another difference is the fact that the wavelengths of 
electric and paramagnetic dipole transitions divide the 
CP potential into two asymptotic regimes: the non- 
retarded regime of distances smaller than these wave- 
lengths and the opposite, retarded regime. The CP po- 
tential follows two distinct l/z A and l/z\ power laws 
in these regimes, respectively. On the other hand, the 
frequency-independent diamagnetic magnetisability lacks 
an intrinsic length scale. As a result, the CP interaction 
with a perfectly reflecting plate follows a retarded l/z A 
power law at all distances. 



IV. TWO-ATOM VAN DER WAALS 
INTERACTION 

Similar to the single-atom case, the body-assisted vdW 
force between two atoms can be derived from the two- 
atom vdW potential U (r A , rg), which is that part of the 
energy shift depending on the positions of both atoms, 



U(r A ,r B ) = AE(r A ,r B ). 



A. Perturbation theory 



(44) 



Again, we assume the atom-field system to be in its 
uncoupled ground-state |0a)|0b)|{0}), but now we have 
to calculate the leading-order two-atom energy shift. We 
begin with two purely diamagnetic atoms, in which case 
the vdW potential follows from the second-order energy 
shift 



A 7 E 



\ H AF 



d 

AF 



Tjd 
BF 



10} 



Eq — Ed 



(45) 



Note that the first order energy shift just leads to the 
sum of two diamagnetic CP potentials as discussed in 
the previous Sect. IIII A[ Due the diamagnetic interac- 
tion Hamiltonian being quadratic in the fields, the vdW 
shift already appears in second order perturbation theory 
rather than in fourth order as for electric and paramag- 
netic dipole transitions. 




|0.4>|0 B >|1,1'> 



FIG. 1: Contributions to the vdW interaction of two dia- 
magnetic atoms. 



The numerator of the term in Eq. (|45[) . when read from 
right to left, represents processes in which the system 
starts from its ground state, goes to a state \<p) due to 
a first atom-field interaction, and finally returns to its 
ground state in the course of a second interaction. As 
the interaction Hamiltonian is quadratic in the magnetic 
field, the intermediate state \<p) must involve the field 
in its ground state or exhibit two field excitations (pho- 
tons). Only the latter case leads to a genuine two-atom 
interaction, 

\4>) = \Q A )\Q B )\l M {r,u),l XH ,(r',u')) (46) 
with the two-photon state being defined by 

|lAi(r, W ),lvi'(r',a;')) - ^/L(r', w')&(r, ^)l{0}>. 

(47) 

The photons must have been emitted by one of the atoms 
and then absorbed by the other. This is schematically il- 
lustrated in Fig. [TJ where the solid lines and the dashed 
lines represent the atoms and the photons, respectively 
(where time progresses in the upwards direction). Note 
that the formal sum in Eq. (|45[) involves sums over A, A', 
i, i' as well as integrals over r, r', oj and lj'. We begin 
with the contribution A2-EV,-) to the energy shift corre- 
sponding to Fig. [Hi), where atom A emits two photons 
and atom B absorbs them, i.e., where the second ma- 
trix element in the perturbative energy shift (|45[) is due 
to the diamagnetic interaction of atom A and the first 
one is due to that of atom B. The required two-photon 
emission matrix element for atom A takes the form 



\H d AF \0) 



-Pa 



Ajk 



x ({0}|/vi'(r',a/)/ Xi (r,a;)4-(r A )4(rA)|{0}), (48) 

where Eqs. ((TTJ) , (|56"]l . and (|47|) have been used. The two- 
photon absorption matrix element {^\H BF \<f)) can be ob- 
tained by taking the complex conjugate of the above and 
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Plate -> 


Perfectly conducting 


Infinitely permeable 


Atom 4- 


Limit — ► 


Retarded 


Nonretarded 


Retarded 


Nonretarded 


Electric 


1 

Z A 


1 

Z A 


I 

Z A 


I 

Z A 


Paramagnetic 


1 
1 

+ - 

Z A 


1 

Z A 


I 

Z A 


1 

Z A 


Diamagnetic 




1 

A 




I 

Z A 



TABLE I: Signs and asymptotic power laws of the ground-state CP potential of an electric, para- or diamagnetic atom with 
a perfectly reflecting plate. 



replacing the labels A by B. After substituting Eq. ([TBI 
for the magnetic field and making use of the commuta- 
tion relations ([2]) and ([3|), the respective term in Eq. (|45] 
reads 



-1 



(o\H d BF \4>)m d AF \Q) 

Eq-E^ 2lj 2 oj' 2 H{lj + lj') 

x tr |/3l • [V A x G* x (r A , r, u) ■ G T x (r B ,r, u) x V B ] 

• pi ■ [V B x G\< (r B ,r', u') ■ G*J(r A , r', w') xV,]}. 

(49) 

At this stage, performing the integrals over r and r' in- 
cluded in Eq. (|4"5l) by means of the integral relation ([22)) 
results in 



27t 2 Jo Jo w + w ' 

xtr [/3^-Im G mm (r j4 ,r B ,a;) ■/3 B -Im G mm (r B , r A) w')] 

(50) 

[Recall Eq. (|27p]. It can easily be seen that the contribu- 
tion (ii) as depicted in Fig. QJii) is exactly the same. We 
hence find 



./ 1 



U dd {r A ,r B ) = - J ^- duj duj 

Jo Jo UJ + UJ' 

x tr [f3 A -lm G mm (r A ,r B ,ui) ■ f3 B -lm G mm (r B , r A , lj')] . 

(51) 

As for the single-atom potential, the result can be sim- 
plified via contour-integral techniques. We first write 



/*°° doj 1 f 
/ — ; — 7 ImG mm (u') = Im / 
Jo w + w Jo 



00 do/ 



w + CJ 



7 G mm (u> ) 



, (52) 

and use the fact that the Green tensor is analytic in the 
upper half of the complex frequency plane including the 
real axis. Hence, we may replace the integral with an 



integral along the positive imaginary frequency axis and 
use the Schwarz reflection principle (|2"0")) to obtain 



du/ 

LJ + Lj' 



huG,„,„iJ) = / d£ 9 i ^ G mm (ig). 

(53) 

Substituting this result into Eq. (|51[) . we next evaluate 
the integral over uj, 



do; • 



lj Im G r 



, c2 



Im 



do; 



wG rom (w) 



(54) 



The integrand in the right hand side has a simple pole 
at uj = i£ in the upper half of the complex frequency 
plane. Again using Cauchy's theorem, we transform this 
integral into a principal value integral along the positive 
imaginary axis, an integral along an infinitesimal half 
circle around the pole, and a (vanishing) integral along 
an infinite quarter circle. The first integral is real, while 
the second integral results in (i7r/2)G OTm (i£), so that 



do; ■ 



lj Im G T 



Im 



V I dv 

o 



vG r 



(iv) 



■ n r 
1 2 

7T _ 



im m (55) 



{V: principal value). After substituting this result to- 
gether with Eq. ([53]) into Eq. ([ST)l. we find the CP po- 
tential of two diamagnetic atoms, 

x tr [(3 A ■ G mm (r A , r B , i£) • f3 B ■ G mm (r B ,r A ,i^)] 
Hf4 



2vr 



x tr [G mm (r^, r B , ' G mm (r B , r a, i£)] (56) 
where the second equality holds for isotropic atoms. 
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(i) 



(ii) 



(iii) 



1 "-v 



B 



B 



B 



FIG. 2: Contributions to the vdW interaction of a diamag- 
netic atom A with an electric atom B. 



Let us consider next the interaction of a diamagnetic 
atom A with an electric atom B. It can be seen easily 
that the leading two-atom energy shift is of third order 
in this case, 



(E^ — E )(E ! p — E ) 



with 



H int =Hi F +H e BF =-±B(r A y$i-B(r A ) - i* B -E(r B ). 

(58) 

The relevant intermediate states \4>) and \^p) can easily 
be determined with the help of the diagrams in Fig. [3] 
Let us begin with diagram [2ji) which correspond to the 
the case where atom A emits two photons and atom B 
absorbs them one after another, so the relevant interme- 
diate states are 

\<f>) = |0 A )|0 B >|lA<(r,«),lA't'(r',« / ))» ( 59 ) 
\il>) = \0 A )\n B )\l x » i »(r",u, , % (60) 

In the two-photon emission matrix element, _ff; n t has 
to be replaced by H AF . This leads, as outlined below 
Eq. gg]), to 

(4>\hUo) 

= -^G* x r (r A ,r,ij)x^ A -p d A -V A xG* x ,(r A ,r',ij'), 

(61) 

while in the two other matrix elements, H lnt has to be 
replaced with H BF . This yields 

fi B -E{r B )\<l>) 

= ^[[vf-Gx(r B ,r,uj)\ i 5y X ,,8 l , l n5{r>-r")5{Lu'-Lu") 

(62) 



{Q\-fr B -E{r B m = -[ t i B n ■Gy,(r B ,r",u%„, (63) 

where the expression (fT4")) for the electric field, as well as 
the commutation relations @ and ([3]) have been used. 
Substitution of these matrix elements into Eq. (I5T[) and 
performing the integrals over r, r' , and r" using the in- 
tegral relation (f2"2"j) . we find the contribution of diagram 
(i) to the energy shift to be 

A 3jB = ^ / du} / dw > ™ 

() Wo Jo + 



X 



[dg» ■ Im K T (r A , r B , w) • /3£ • Im K (r A , r B , w' ) • <] . 



(64) 



where 



K(r,r',w) = VxG(r,r» 



(65) 

The contribution to the energy shift from diagrams (ii) 
and (iii) result in terms similar to Eq. (|64j) except that 
the energy denominator (oj +lu')(u b + u) has to be re- 
placed with (—uj b — uj')(lj b + lo) for diagram (ii) and 
with {uj + ui')(tjJ B + lu') for diagram (iii). Summing all 
contributions, wc find 



A,£ = 



2//; 



da; 



* 2 Jo Jo (wS+w)(w + w')(wS+a;') 
x [dg» ■ Im K T (r A , r B , w) • (3 d A ■ Im K (r A , r B , w' ) • <] . 

(66) 

By transforming the u/ integral by means of contour- 
integral to run along the positive imaginary axis, 



, u'ImK(w') 

o — ( w + w')K + w) 



do; 



■ ?(a,g + M )K(iQ 
o d V + a [ (wS )2 + ^]' lb ° 

and then performing the u integral in the way explained 
above Eq. (fool) , one obtains 

U de (r A ,r B ) = ^J™dte 

x tr [/3 A ■ K [r A , r B , i£) • ol b (i£) • K T (r A , r B , if)] 

t 2 /-oc 

= dH 2 P A a B m 
2tt Jo 

x tr[K(r A , r B , iO • K T (r A , r B , , (68) 

where the definition AMI) has been used and the second 
equality holds for isotropic atoms. Our result (|68)l for the 
vdW potential of a diamagnetic atom ^4 with an electric 
atom B may be written in the form (/3 A = c 2 ol a ) 



U de (r A ,r B ) = -- — j / d£ 
2tt£o Jo 



tr[a A -G me (r A ,r B A£)-a e B (i£)-G em (r B ,r A ,i£)] 



fc p CO 



tr [G me (r A , r B ,i£)-G em (r B ,r A , i£)] , (69) 
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where 



G me (r,r» = — K(r,r',w) - — VxG(r,r», (70) 

c c 

G em (r,r',u;) = ^K T (r>,w) = -G(r,r»x$'. 

c c 



(71) 



Obviously, the interaction between an electric atom A 
with a diamagnetic atom B is given by the right hand side 
of Eq. (|68| with the labels A and B being interchanged. 

The vdW potential of a diamagnetic atom A with a 
paramagnetic one B can be calculated in a completely 
analogous way. The lowest-order energy-shift is again 
given by third-order perturbation theory with the pos- 
sible intermediate states given again as depicted in the 
diagrams in Fig. [3J The result is (J3 A = c 2 a A ) 



h f°° 
U dP (rA,r B ) = -^J o df 



x tr [a A ■ G mm {r A , r B , if) • of B (if) • G mm (r B , r a , if)] 

-n 



2^1 



tr [G mm (r A , r B , if ) -G mm (r B , r A , if)] . 



(72) 



We have thus generalised the body-assisted vdW po- 
tential to allow for atoms with a diamagnetic response. 
Recalling the previously derived interaction of two para- 
magnetic atoms [Hj], we again observe that the gen- 
eralization can be achieved by replacing the paramag- 
netisability of each atom with the total magnetisability, 
(3p(lo) h-> (3(w) = f3 p (oj) + (3 d . The complete vdW po- 
tential of two atoms with nontrivial electric, para- and 
diamagnetic properties can thus be given as 



U(r A ,r B ) = 



X.X' —e.m 



Uxx>(r A ,r B ), (73) 



with 

U\\>(r A ,r B ) 



2tt£o Jo 



df 



tr [ a A (if) ' G XX' ( r A , r B , if) • a.% (if) • G y x (r B , r A , if)] 



2ne 2 



dfal(if)a B (if) 



o Jo 



tr[G X x' (r A , r B , ig)-Gx\(r B , r A , if)] , 



(74) 



Again, this implies that previous results for local-field 
corrected vdW potentials [18j remain valid for diamag- 
netic atoms. 



B. Application: Two atoms in free space 

As the simplest example for the two-atom interaction, 
let us consider two isotropic atoms A and B interacting 



with each other in free space. The Green tensor reads 
(see e.g. Ref. E3) 

G<°>(rA,r B ,iO - ^3 [/af/c)l- 5 af/c)e ie; ]e-'«/ c 

(75) 

(r A y^r B ), where l = r A r B , I — \l\, ei = l/l, and 

f(x) = 1 + x + x 2 , g(x) = 3 + 3x + x 2 . (76) 

Evaluating the curls as contained in definitions (f2~7)) . ([70)1 
and (|7ip . one easily finds 

GW(r A ,r B ,iO = G%l(r A ,r B ,iO 

_ |__(o) 



G^(r A ,r B ,if), (77) 



G^l (TA,r B , if) = -GW(r A ,r fl , if) 



:(l + /f/c)e^«/ c e ; xl. (78) 



Substituting these into Eqs. (174")) and (175]) we obtain a 
total free-space vdW potential 



U(r A ,r B ) = U(l) 



16tt 3 



x i 76" / d ? [c 4 aA(if)a B (if)+/3A(if)/3 B (if)]/ i i(?f/c) 



+ 7J / dff 2 [a A (if)^ B (if)+/3^(if)a B (if)]/i 2 (Zf/c) 



7 4 / 

f JO 



where 



hi(x) = (3 + 6.t + 5a: 2 + 2x 3 + x 4 )e~ 2x 
h 2 {x) = (1 + 2x + x 2 )e~ 



(79) 



(80) 



As expected from the results of Sect. IIV A[ this result 
has the same form as the previously deriv ed p otential 
between electric and paramagnetic atoms [la ], except 
that here the total magnetisability appears in place of 
the paramagnetic one. 

It is of particular interest to inspect the behaviour 
of the interaction potential in the nonretarded/retarded 
limits, where the atom-atom separation is small/large 
compared to the respective atomic wavelengths. We 
focus here on the contribution of the atomic diamag- 
netisabilities to the potential, given by Eq. (179")) with 
/3(if) -> (3 d . 

The interaction potential between a diamagnetic atom 
A and an electric atom B, using the explicit expression 
for the polarisability (|3~4l in Eq. (f79")> . is found to be an 
attractive potential in the form 



U de {l) 



^^cMf/ df , ( 

o f^ 



24tt 3 ; 4 ^ 

k 



e 



^(/f/c). 



(82) 
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To achieve the limiting cases mentioned above, we note 
that in the ^-integral, ujg <C £,{^b ^> £) holds in the non- 
retarded (retarded) limit. This leads to a Z~ 5 -dependent 
and a Z _7 -dependent potential for the nonretarded and 
retarded limits, respectively, as follows 



96tt 3 / 5 



Ofc |2 



(83) 



u T de {i) 



96tt 3 Z 



1Mb 



Ofc 1 2 



|/3l|a B (0). 



The diamagnetic-paramagnetic part of the two-atom 
interaction in free space is seen to be repulsive. It reads 



24tt 3 / 6 



E>£K»T/ &H—— 1 hM/c), 

k Jo e+u% 

(85) 

where Eq. (1351) is used for the paramagnetisability of 
atom £?. For the nonretarded limit Eq. ([85"T) exhibits a 
Z~ 6 -dependence, 



16tt 2 Z 6 



E 



m 



fe |2 



16ir 2 l e 



(86) 



while in the retarded limit it tends to a I 7 -dependent 
potential, 



uui) 



96tt 3 Z 7 



E 



m 



Ofc 12 



(87) 

Finally, the diamagnetic-diamagnetic two-atom in- 
teraction in free space shows a unique attractive in - 
dependence for any arbitrary range of atom-atom sep- 
aration, due to the frequency- independence of diamag- 
netisabilities. As can be obtained from Eq. ([751) . it is 
given by 

In Table m we compare the signs and asymptotic 
power law of the vdW potentials involving diamagnetic 
atoms with the known results for purely electric or para- 
magnetic atoms [l8[. We observe that the replacement 
of a paramagnetic atom with a diamagnetic one leads 
to a sign change of the vdW interaction as predicted by 
the Lenz rule. In addition, the frequency-independence 
of the diamagnetic magnetisability leads to new asymp- 
totic power laws, e.g., l/l 5 for the nonretarded electric- 
diamagnetic interaction. This result is in between the 
l/l 6 and l/l 4 asymptotes found for the electric-electric 
and electric-paramagnetic cases. 



C. Comparison with Microscopic QED 

It is instructive to compare the results obtained in this 
section for the vdW dispersion interaction between atoms 



in free space when either one or both species is diamag- 
netic using body-assisted fields, with that derived using 
microscopic QED. In this second approach it is com- 
mon to use the well-established molecular QED theory 
[I5I [l6j . For two mutually interacting particles A and 
B in vacuum coupled to the radiation field, the total 
Hamiltonian is written in the form of Eq. (|13[) with 
given by Eq. ([1} . In the multipolar coupling scheme, the 
Hamiltonian operator for the free electromagnetic field is 
expressed as 



Hp 



d 3 7 



1 



— d 2 {r) + —b 2 (r) 

£0 Ho 

= E[a (A)t W« (A) (fc) 

k,\ 



hw (89) 



where d(r) and b(r) are second quantized microscopic 
electric displacement and magnetic field operators, re- 
spectively. They are commonly written as a mode sum 
in terms of vacuum boson annihilation and creation oper- 
ators a^(k) and a^"(k) for a photon of wave vector k, 
polarization index A, and circular frequency oj = ck 7 as 
in the second equality of Eq. (|59")l . Analogous to Eq. © 
with Eqs. (HI)-®, the atom-field coupling Hamiltonian 
for the dispersion interactions of interest is 



h sf = - £ o l fre d ( r s) - rh ?- fe ( r e) - 



\ - T Q 



8m c 



[f Q x6(r e )] 2 

It is worth noting the explicit appearance in the last 
two expressions of the electric displacement field oper- 
ator. This is a direct consequence of adopting the mul- 
tipolar framework, where the field momentum canoni- 
cally conjugate to the vector potential is proportional 
to d(r) instead of to the electric field itself. This is a 
common feature of the microscopic and macrocopic ap- 
proaches. In the latter, the quantised field E(r) also 
refers to the Power-Zienau transformed electric field op- 
erator that has to be regarded as displacement field with 
respect to the atomic polarisation [12j ]. 

Dispersion potentials between a diamagnetic atom and 
an electrically polarizable and a magnetically suscepti- 
ble atom, and between two diamagnetic atoms may be 
computed in a manner similar to that already detailed. 
As before, the last interaction is obtained via second- 
order perturbation theory and is depicted in Fig. [TJ while 
the former two occur in third-order and are described by 
Fig. [21 Matrix elements are evaluated using unperturbed 
product atom-field states \n^)\m(k, A)) = \n^;m(k, A)) 
where a number state representation is used to signify 
the number of photons present, in this case m. 

As expected, results identical to Eqs. (|52|) . ([55} and 
([55)1 are found for the vdW dispersion potentials between 
a diamagnetic and an electric atom, a diamagnetic and a 
paramagnetic one, and between two diamagnetic atoms, 
valid for the entire range of separation distance vector 
I beyond wave function overlap of the two centres and 
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Atom A -t 


Electric 


Parama 


gnetic 


Diama; 


metic 


Atom B 1 


Limit — > 


Retarded 


Nonret. 


Retarded 


Nonret. 


Retarded 


Nonret. 


Electric 


1 

~F 


1 

~¥ 


1 

-I 


1 

A 

+ l± 


1 

~¥ 


1 

~F 


Parama 


gnetic 


1 

+ ¥ 


1 

+ ¥ 


1 

~F 
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~¥ 


l 

+ ¥ 


l 

+ ¥ 


Diama; 


metic 


i 


l 


l 

+F 


i 

+ ¥ 


~i 


L 

7 



TABLE II: Signs and asymptotic power laws of ground-state vdW potentials of electric, para- or diamagnetic atoms in free 
space. 



extending out to infinity. Explicit details of this cal- 
culation may be found in Ref. [25j]> and the respective 
energy shifts are given by Eqs. (3.7), (3.20) and (3.43) 
of that paper. Furthermore, identical near- and far-zone 
asymptotic limits clearly follow in each case. The connec- 
tion with the macroscopic approach can be established 
by computing all necessary field expectation values and 
evaluating them using the free-space Green tensor ([75]). 

V. DISCUSSION AND SUMMARY 

We have calculated CP and vdW potentials of atoms 
with nontrivial diamagnetic properties in the presence of 
arbitrary magnetoelectric bodies on the basis of macro- 
scopic QED and leading-order perturbation theory. The 
nonlinear interaction generating the diamagnetic interac- 
tion is quite different from the paramagnetic coupling and 
it leads to different perturbative orders. Nevertheless, 
we have found that diamagnetic atomic properties lead 
to dispersion potentials which formally resemble those 



of paramagnetic atoms where the diamagnetic magneti- 
sability appears in place of the paramagnetic one. We 
have explicitly shown this correspondence for one- and 
two-atom potentials, but it is expected to hold for multi- 
atom vdW potentials as well. 

However, the fact that the diamagnetic magnetisabil- 
ity is negative and frequency-independent leads to dia- 
magnetic potentials that differ in signs and power laws 
from their paramagnetic counterparts. Diamagnetic dis- 
persion interactions carry the same sign as the well- 
known electric potentials, which implies that diamag- 
netism alone cannot be used to realise repulsive po- 
tentials. The unique power laws resulting from the 
frequency-independence imply that diamagnetic poten- 
tials have their strongest influence at short range. 
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